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L  Background 


Electrons  produced  by  various  ionization  processes  can  play  an  important  role  in 
detennining  spectral  radiative  intensities  and  odier  characteristics  of  hardbody  flowfields  and 
rodket  fdumes.  Robust  modelling  of  the  radiaticm  emitted  from  such  flows,  as  well  as  of  other 
idiysical  and  chemical  pn^ierties,  must  include  collisicmal  processes  involving  these  low-eneigy 
electrais  with  the  molecular  species  inesmit  Although  the  collision  energies  are  low,  non- 
equililnium  behaviOT  in  hardbody  flowfields  can  kad  to  overshoots  resulting  in  temperatures  in 
die  12,000  K  range.  Bectrons  in  the  tail  of  these  oiergy  distributions  can,  and  are  believed  to, 
cause  dectronic  excitation  of  species  in  these  flowfields.  Cross  sections  for  such  threshold 
excitation  by  electrons  with  just  sufficient  energy  to  drive  the  process  are  generally  inaccessible 
experimentally  and  the  data  base  of  these  cr<»s  sections  is  essentially  nonexistent  Furthermore, 
impact  oi  very  low-mm'gy  dectrons,  present  in  relatively  high  number  densities  in  these  low- 
temperature  plasmas,  on  electronically  exdted  metastable  states  can  lead  to  significant 
dqx^Milation  of  these  states.  Very  little  is  known  about  such  cross  sections  for  electron  impact 
on  metastable  states  other  than  they  are  expected  to  be  quite  large. 

In  summary,  many  of  the  cross  sections  for  electron-molecule  collisions  needed  for . 
robust  modelling  of  the  properties  of  hot  interacting  flow  fields  associated  with  hardbodies  and 
idnmes,  such  as  radiative  signatures,  are  not  known.  Use  of  unreliable  or  guessed  values  of  these 
cross  sections  can  lead  to  incorrect  models  of  the  underlying  processes  responsible  for  radiative 
agrutures  and  their  behavior  with  velodty  and  altitude.  Models  based  on  such  urueliable  cross 
aectkms  can  fail  seriously  outside  of  the  limited  range  where  some  validation  may  have  been 
qrparenL 

One  strategy  for  obtaining  these  cross  sections,  particularly  those  which  ate 
eiqrerimentally  iruiccessible,  is  to  calculate  thetiL  The  major  hindrance  to  such  a  strategy, 
however,  has  been  the  huge  computational  demands  associated  with  the  solution  of  the  equations 
governing  electton-molecule  collisions  at  low  impact  enogies.  In  fact,  it  is  widely  recognized 
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thtt  die  coet-to-peffonnanoe  ratio  of  convendooal  vector  supercompiiteis  such  as  the  CRAY  Y- 
MP  has  made  such  an  ^^ifoadi  impracticaL 

n.  Technical  Summary 

A.  Objective 

The  objective  of  this  efimt  was  to  exploit  a  thecwedcal  fonnulation  and  computatitmal 
methodology  which  we  have  ^wcifically  devel(^)ed  for  such  applkadons  to  calculate  die  cross 
sections  few  dectnxi  impact  mtdtadon  of  mc^ules  such  as  N2, 02.  NO,  OH,  CN,  and  H2O. 
These  cross  sections  are  needed  for  robust  modelling  of  the  ultraviolet  emission  occurring  in  die 
bow  diode  of  a  ballistic  missile.  Such  ultraviolet  signatures  have  a  clear  potmitial  for  detection 
and  tracking  of  ballistic  and  theater  missiles.  The  calculation  of  these  cross  sections  for 
molecules  is  a  computationally  intensive  quantum-mechanical  problem.  These  cranputaticmal 
demands  and  complexities  have  sevmely  limited  {sogiess  in  such  studies  to  date.  As  we  will  see 
diordy,  one  cf  our  key  accomplishments  has  been  our  ability  to  successfully  exploit  the  high- 
petfrnmance  and  cost-effective  computing  provided  by  massively  parallel  supeiccnnputms  to 
meet  die  computational  demands  of  these  calculatimis.  These  parallel  supercomputers  achieve 
very  high  aggn^ate  qieeds  and  large  memtny  by  harnessing  die  power  of  cmnmercially 
available  microprocessois  assembled  in  a  scalable  architecture.  Such  parallel  ctHnputers  can  be 
expected  to  revolutionue  our  ability  to  model  complex  physical  and  chemical  systems. 

B.  Summary  of  Accomplishment 
(a)  Background 

To  put  our  accomplishments  in  perspective  it  is  essmitial  to  review  a  few  salient  features 
of  our  fonmilation  of  the  electron-molecule  collision  problem  and  the  computational 
methodology  used  in  diese  studies.  As  in  many  complex  physical  systmns,  we  do  not  attempt  to 
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die  relevant  wave  eqn^on  governing  dieae  «dlisioos  bat,  instead,  employ  a  variatkmal 
princ^de  wUdi  reduces  the  problem  to  one  of  stdving  ^sterns  cd’ linear  equations  for  the 
quantity  of  direct  idiysical  inieiesL  In  these  a|q;>lications  diis  quantity  is  the  scattering  amplitude 
fiom  which  die  anoss  section  is  determined.  We  use  a  multichannel  extension  of  the  variational 
prine4de  for  die  scattering  amplitude  originally  introduced  by  Schwinger.^  This  Schwinger 
multichannel  variational  procedure  was  specifically  formulated  for  applicaticms  to  electron- 
molecale  ooUiakms.^  A  key  feature  of  this  multichannel  variational  procedure  is  diat,  as  in  the 
original  Schwinger  principle,  die  trial  wave  function  need  not  satisfy  scattering  boundary 
conditums;  in  particular,  square-integrable  functicms  such  as  Cartesian  Gaussians,  commonly 
used  in  molecular  electronic  structure  calculations,  may  be  employed. 

Use  cf  diis  variaticmal  principle  leads  to  a  systmn  of  linear  equations 

Ax  =  b.  (1) 

vrimse  solutitms  x  yield  the  scattering  amplitude.  In  Eq.  (1)  A  is  a  complex  symmetric  matrix 
and  b  and  x  are  rectangular  complex  matrices.  The  dements  of  the  matrices  A  and  b  have  the 
form 

a,=<^i(i^-p]h+vp-vg<;>vi0,>  (2) 

and 

iv.=<^ilVI«».(r„r„-r,)exp(ik.-r,J>.  o) 

In  Eq.  (2)  A  and  dj  are  (N  -f  l)-electron  Slater  determinants  in  which  the  trial  scattering  wave 
functum  for  the  composite  system  of  electron  (rius  molecule  is  expanded,  H=E  -  H,  wfame  H  is 
die  Hamiltonian  for  die  incident  electron  plus  motecular  target  with  N  electrons,  V  is  the 
potential  energy  operatm’  ftx  die  interactimi  of  die  incident  electron  with  die  target,  P  is  a 
projection  operator  vriiich  selects  opon  channels  (energetically  allowed  outcomes  of  the 
ctdliaion),  and  0]^^  is  die  intmacticm-fm  (3ieen's  function  with  tqipropriate  boundary 
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conditioas,  pcojected  onto  the  (^>en  channds.  In  Eq.  (3)  is  an  dgastate  of  the  molecular 
tttset  with  energy  E.  andmcp  (ika>rN^i)descrU)es  a  plane  wave  pn^gating  in  the  direction 
given  by  ka.  The  angle  toidc^  in  Eqs.  (2)  and  (3)  indicate  integration  over  the  coordinates  of 
aD  electrons. 

Construction  of  die  matrices  A  and  b  requires  die  evaluation  of  several  different  types  of 
im^prals.  The  difficulty  w  ease  of  evaluadonttf  die  matrix  elements  arising  in  a  vaiiadtHial 
piincqile  is  cleariy  a  maj<M’  factor  in  die  choice  of  mediod.  In  our  variational  principle  the  trial 
wave  function,  Le..  die  ^  and  of  Eqs.  (2)  and  (3),  need  not  satisfy  scattering  boundary 

ctmditicms.  In  fact,  if  we  exploit  this  flexibility  and  expand  the  molecular  orbitals  in  the  Slato- 
d^erminants  and  in  Cartesian  Gaussian  functions  of  the  form 

o(r;/,in,n,R,{:)=N.(x-X)'(y-y)“  (z-Z)*exp(-flr-RPX  (4) 

where  N«  is  a  nonnalizatitm  omstant  and  R  s  (X,Y,Z),  all  integrals  associated  with  the 
Hamiltonian  (H)  and  potential  (V)  terms  in  Eqs.  (2)  and  (3)  can  be  evaluated  analytically. 
Effident  techniques  and  algorithms  have  been  developed  fOT  their  evaluation.  However,  the 
integrals  involving  VGjT^  V  have  no  known  arulytic  form  and  ate  by  far  the  most  difficult  and 
computationally  intmisive.  Quadrature  of  smne  fcmn  must  be  employed  to  evaluate  these  trains. 

(b)  Computational  Breakthrough 

The  occurrmce  of  terms  containing  the  Gierai's  function  in  variational  principles  based 
<m  integral  equations  has  oftrai  discouraged  their  use,  in  ^te  of  known  braiefits  of  these  methods 
such  as  die  flexilnlity  of  using  trial  fonctimis  udiich  do  not  satisfy  scattering  boundary  conditions 
and  dteb  favorable  convergraice.^  A  signifk^t  accranplishment  of  our  wrak  has  been  the 
develoimirait  of  a  strategy  fra  evaluation  ci  diese  Green's  function  natiix  elranraits  on  parallel 
computrars  consisting  of  huncbeds  of  commercially  available  microprocessors.  This  strategy  has 
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been  uccesriuHy  implNneated  (m  the  Intel  Toadist(»e  Delta  S]rstem.  a  multictnuputn’ 
containing  S12  Intel  j860micn^»ooessocs.  Runs  to  date  (» the  Intel  Delta  show  sustained 
perfonnanoe  oi  aevenl  billions  of  arithmetic  operadcms  per  sectmd  (Gigaflops)  in  calculatimu  of 
the  cross  aectkms  of  interest  to  us.  This  petfonnance  far  exceeds  that  seen  (mconventumal 
siqieicoiiyutera  such  as  the  CRAY  Y-MP  by  mcxe  than  a  factor  of  1(X)  and,  with  increasing 
micrc^irooessor  qweds  and  scaling  of  these  systems,  can  be  expected  to  increase  significantly  in 
die  near  future.  This  deveU^mtent  opens  up  endrefy  new  vistas  in  ampuuaional  studies  of  Ais 
kind  and,  m  fact,  is  an  earty  indication  of  Ae  revolutionary  impact  duit  scalable  high- 
performance  com/mting  can  be  expected  to  have  on  our  alnUty  to  model  com/dex  physical 
systems. 

hi  our  quadrature  of  diese  matrix  elements  involving  VGp^W,  the  i»incipal  stq>  is  the 
evaluation  and  subsequent  transformation  of  a  large  number,  e.g.,  10^  or  10*^  of  two-dectrtm 
integrals  of  die  form 

J  d’rjj  d%a(r,)h(r,)-^y(r,)exp(ikr,), 

^12 

involving  three  Cartesian  Gaussian  functions  (a,^,  and  y)  and  a  plane  wave  exp  QkT). 

AldKmgh  diese  integrals  can  be  evaluated  analytically  via  a  few  diousand  lines  of  Fortran,  their 
evaluation  can  be  computadtHiaUy  demanding  because  the  numb^  of  such  integrals  is  large. 

These  two  cmnputationally  intensive  tasks  (0  the  evaluation  of  a  very  large  number  of 
dementtty  integrals  describing  the  interactitm  between  etectrons,  and  (I’O  die  transformaticHi  of 
diese  elementary  integrals  mto  the  many-^tectrcm  integrals  appearing  in  die  variadtuud 
eqHesskm,  are  ideally  suited  for  impkmentatitm  on  distributed-memory  paralld  ctunputers. 
These  machines  consists  of  fairly  pownful  cmnmercially  available  microiMrocess(»s  with 
manory  which  <^)erate  autcMKnnously  and  ctmununicate  with  each  other  via  a 
mrssagfr-pasiing  network.  The  evaluatkm  of  the  large  sets  of  two-electron  integrals  is  readily 
paraPdiaed:  one  simply  has  eadi  micit^nocessor  independendy  calculate  a  different  subset  of 
integrals.  Such  a  parallel  strategy  is  especially  appealing  since  die  integral-evaluation  algOTithms 
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involve  recorsive  steps  that  make  vectorizatiiMi  on  scquentiil  supercomputers  difficult  The 
tnmsfonnation  of  diese  integrals  cm  be  organized  as  a  series  of  large  distributed  matrix 
nwihiplications  between  the  global  array  of  integrals  stored  in  the  microprocessors  and  a 
distributed  transfomation  matrix  that  is  easy  to  construct  Since  thoe  ate  gmierally  too  many 
integrals  to  hcdd  in  memory  simultaneously,  we  treat  y  and  the  magnitude  of  k  in  (S)  as 
sequential  indices.  The  relative  importance  of  these  two  tasks  depends  (m  the  problem  at  hand, 
but  poformance  in  each  step  far  exceeds  diat  achieved  by  the  original,  sequential  program.  A 
large  productitm  tun  on  the  full  Intel  Delta  can  sustain  througlq>ut  of  about  2J  GFLOPS, 
apprmdmately  100  times  tiu  performance  of  tite  original  sequential  program  on  the 
CRAYY-MP. 


(c)  Most  Significant  Results 

We  have  exploited  this  parallel  computational  procedure  and  processor  time  on  the  Intel 
Ddta  System  to  carry  out  smne  of  the  first  robust  calculaticms  of  the  cross  sections  for  excitation 
of  die  low-lying  excited  states  of  N2,  CO,  H2O,  and  other  molecules  by  electrons  with  impact 
energies  close  to  threshold.  Although  this  threshold  excitation  region  is  most  relevant  in  the 
modelling  of  these  flowfields,  it  is  also  generally  not  accessible  expoimentally. 

Our  calculated  cross  sections  for  excitation  of  the  A  state,  the  lowest-lying 
metastable  state  of  N2  with  a  lifetime  of  ~1  show  a  rapid  onset  at  threshold  and  agree  well 
with  the  single  set  of  measured  values  which  ate  available  at  a  few  electron  volts  above 
ditediold.^  No  measured  values  ate  available  closer  to  threshold.  Due  to  its  very  long  lifetime 
tins  metastable  state  plays  an  important  role  in  modelling  the  radiative  signatures  from 
flowfields.  Fbrccmvenience,  these  cross  sections  are  shown  in  the  attached  figure.  Our 
calculated  values  provide  the  most  reliable  estimates  of  the  cross  sections  for  excitation  of  the 
A  state  orN2  by  electnm  impact  at  near-threshold  energies.  Cross  sections  for  electron 
inqwct  excitation  of  many  otho'  low-lying  excited  states  of  N2  have  also  beoi  calculated. 
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Agieaiient  with  measmed  values  of  these  ctoss  sections  is  quite  encouraging  in  ^ite  of  the 
significant  inconsistencies  anumg  the  different  data  sets. 

Similar  cakulaticMis  have  been  carried  out  for  the  cross  sections  for  electnm  impact 
excitation  of  die  lowest  lying  m^astable  a  ’ll  state  of  CO  which  is  the  upper  levd  of  the 
Cameron  bands  of  CO.^  The  cross  sectitms  also  ^ow  a  very  nqiid  onset  at  threshold  in 
agreement  with  values  obtained  from  tr^iped  electron  measuremoits  of  thirty  years  ago.^  Due  to 
well-known  difiBculties  associated  with  the  technique,  trapped  electron  measurements  were  not 
genenlly  viewed  as  reliable.  Cross  sections  for  electron  impact  excitation  of  many  odier  low- 
lying  excited  sbues  of  CO  have  also  been  calculated.^  These  calculations  are  being  extended  to 
obtain  the  cross  sections  fOT  excitation  of  the  a ’ll  ^ate  of  CO  via  electron  impact  on  high 
vitoadonal  levds  of  die  ground  state.  These  cross  sections  could  play  a  significant  role  in 
modelling  systons  with  substantial  peculations  of  molecules  in  high  vibrational  levels. 

We  have  also  cmnpleted  calculations  of  the  cross  sections  fw  excitation  of  the  B  ’ll, 

state  of  Nj  by  electron  impact  on  the  metastable  A  ’x^  state.  Electron  impact  on  metastable 
stales  can  lead  to  significant  depopuladem  of  excited  electronic  states  in  die  low-temperature 
plasmas  occurring  in  flowfields.  The  small  energy  gic^  between  excited  electronic  states,  the 
high  density  of  such  excited  states,  and  the  larger  number  densities  of  electrons  with  sufficient 
energy  to  drive  such  excitations  all  contribute  to  the  importance  of  this  mechanism.  To  our 
knowledge,  these  calculaticms  provide  the  only  available  estimates  of  cross  sections  for  electron 
impact  excitation  of  metastable  molecular  nitrogen.  Further  studies  along  these  lines  are  under 
way  for  N2  and  CO. 

We  have  also  exploited  our  fmmulation  and  computational  procedure  to  obtain  cross 
sections  for  exdtatitni  and  dissociadem  of  H2O  by  electron  impact,  i.e.,  e  +  H2O H-i- OH.  At 
die  time  oi  publication,  these  were  the  first  results  of  ab  inido  calculadons  of  such  cross  secdons 
for  a  ptdyafomic  mokcule.  Estimates  of  the  cross  secdon  for  electron  impact  dissociadon  of 
water  based  on  die  results  of  these  calculadons  diffo^  dramadcally  from  the  recommended  values 
availalde  ki  die  litmature.'^ 
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(d)  Siwunary  and  Outlook 


The  high-perf(»maiice,  cost-effective  computing  provided  by  massively  parallel 
computers  can  be  expected  to  have  a  significant  impact  on  our  ability  to  simulate  complex 
(diysical  and  engineering  systems.  Such  capabilities  will  strengthen  the  competitiveness  of  U.S. 
defoise  and  cmnmercial  industries.  However,  robust  simulations  of  systems  such  as  flowfields 
over  hanfiK)dies  and  plumes  must  increasingly  rely  on  an  adequate  data  base  for  the  fundamental 
chemical  and  fdiysical  fnocesses  occurring  in  these  systems.  Much  of  this  data  base  is  dtho’ 
unavailable  (V  experimentally  inaccessible. 

Our  use  of  parallel  ccmiputers  in  diese  studies  of  electron  collisions  with  molecules  for 
modelling  of  flowfields  is  an  eariy  example  of  the  significant  role  that  scalable  high-performance 
C(»nputing  can  play  in  meeting  such  data-base  needs  for  several  defense  and  dual-use 
^tplications.  With  the  projected  scaling  of  these  machines  to  much  larger  numbers  of 
micrt^HOcessors  and  improvements  in  hardware  and  algorithms,  performance  levels  weU  in 
excess  of  diat  achieved  (m  the  eariy  versions  of  the  machines  used  in  these  studies  can  be 
expected.  In  summary,  these  machines  will  certainly  have  a  revolutionary  impact  on  our  ability 
to  model  cmnplex  systems. 
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